A Search for Time-Dependent B® — B® Oscillations 

Using Exclusively Reconstructed Df Mesons 

Kenji Abe/ 15 ) Koya Abe/ 24 ) T. Abe/ 21 ) I. Adam/ 21 ) H. Akimoto/ 21 ) D. Aston/ 21 ) 
K.G. Baird/ 11 ) C. Baltay/ 30) H.R. Band/ 29) T.L. Barklow/ 21 ) J.M. Bauer/ 12 ) 

G. Bellodi/ 17 ) R. Berger/ 21 ) G. Blaylock/ 11 ) J.R. Bogart/ 21 ) G.R. Bower/ 21 ) J.E. Brau/ 16 ) 
M. Breidenbach/ 21 ) W.M. Bugg/ 23 ) T.H. Burnett/ 28 ) P.N. Burrows/ 17 ) A. Calcaterra/ 8 ) 

R. Cassell/ 21 ) A. Chou/ 21 ) H.O. Cohn/ 23 ) J.A. Coller/ 4 ) M.R. Convery/ 21 ) R.F. Cowan/ 13 ) 

G. Crawford/ 21 ) C.J.S. Damerell/ 19 ) M. Daoudi/ 21 ) N. de Groot/ 2 ) R. de Sangro/ 8 ) 
D.N. Dong/ 21 ) M. Doser/ 21 ) R. Dubois/ 21 ) I. Erofeeva/ 14 ) V. Eschenburg/ 12 ) S. Fahey/ 5 ) 
D. Falciai/ 8 ) J.R Fernandez/ 26 ) K. Flood/ 11 ) R. Frey/ 16 ) E.L. Hart/ 23 ) K. Hasuko/ 24 ) 
S.S. Hertzbach/ 11 ) M.E. Huffer/ 21 ) M. Iwasaki/ 16 ) D.J. Jackson/ 19 ) P. Jacques/ 20 ) 
J.A. Jaros/ 21 ) Z.Y. Jiang/ 21 ) A.S. Johnson/ 21 ) J.R. Johnson/ 29 ) R. Kajikawa/ 15 ) 
M. Kalelkar/ 20 ) H.J. Kang/ 20 ) R.R. Kofler/ 11 ) R.S. Kroeger/ 12 ) M. Langston/ 16 ) 
D.W.G. Leith/ 21 ) V. Lia/ 13 ) C. Lin/ 11 ) G. Mancinelli/ 20 ) S. Manly/ 30 ) G. Mantovani/ 18 ) 
T.W. Markiewicz/ 21 ) T. Maruyama/ 21 ) A.K. McKemey/ 3 ) R. Messner/ 21 ) K.C. Moffeit/ 21 ) 
T.B. Moore/ 11 ) M. Morii/ 21 ) D. Muller/ 21 ) V. Murzin/ 14 ) S. Narita/ 24 ) U. Nauenberg/ 5 ) 

H. Neal/ 30 ) G. Nesom/ 17 ) N. Oishi/ 15 ) D. Onoprienko/ 23 ) R.S. Panvini/ 27 ) C.H. Park/ 22 ) 

I. Peruzzi/ 8 ) M. Piccolo/ 8 ) L. Piemontese/ 7 ) R.J. Piano/ 20 ) R. Prepost/ 29 ) 
C.Y. Prescott/ 21 ) B.N. Ratcliff/ 21 ) J. Reidy/ 12 ) PL. Reinertsen/ 26 ) L.S. Rochester/ 21 ) 
PC. Rowson/ 21 ) J.J. Russell/ 21 ) O.H. Saxton/ 21 ) T. Schalk/ 26 ) B.A. Schumm/ 26 ) 
J. Schwiening/ 21 ) V.V. Serbo/ 21 ) G. Shapiro/ 10 ) N.B. Sinev/ 16 ) J.A. Snyder/ 30 ) 
g^; H. Staengle/ 11 ) A. Stahl/ 21 ) P. Stamer/ 20 ) H. Steiner/ 10 ) D. Su/ 21 ) F. Suekane/ 24 ) 

A. Sugiyama/ 15 ) S. Suzuki/ 15 ) M. Swartz/ 9 ) F.E. Taylor/ 13 ) J. Thorn/ 21 ) E. Torrence/ 13 ) 
T. Usher/ 21 ) J. Va'vra/ 21 ) R. Verdier/ 13 ) D.L. Wagner/ 5 ) A.P Waite/ 21 ) S. Walston/ 16 ) 
A.W. Weidemann/ 23 ) E.R. Weiss/ 28 ) J.S. Whitaker/ 4 ) S.H. Williams/ 21 ) S. Willocq/ 11 ) 
R.J. Wilson/ 6 ) W.J. Wisniewski/ 21 ) J.L. Wittlin/ 11 ) M. Woods/ 21 ) T.R. Wright/ 29 ) 
R.K. Yamamoto/ 13 ) J. Yashima/ 24 ) S.J. Yellin/ 25 ) C.C. Young/ 21 ) H. Yuta/ 1 ) 

(The SLD Collaboration) 



< - 1 ) Aomori University, Aomori , 030 Japan, 
( 2 ) University of Bristol, Bristol, United Kingdom, 
^Brunei University, Uxbridge, Middlesex, UB8 3PH United Kingdom, 
^Boston University, Boston, Massachusetts 02215, 
( 5 ) University of Colorado, Boulder, Colorado 80309, 
^Colorado State University, Ft. Collins, Colorado 80523, 



1 



^INFN Sezione di Ferrara and Universita di Ferrara, 1-44100 Ferrara, Italy, 
&INFN Lab. Nazionali di Frascati, 1-00044 Frascati, Italy, 
Johns Hopkins University, Baltimore, Maryland 21218-2686, 
^Lawrence Berkeley Laboratory, University of California, Berkeley, California 94720, 
( n ) University of Massachusetts, Amherst, Massachusetts 01003, 
( 12 ) University of Mississippi, University, Mississippi 38677, 
^Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, 
^Institute of Nuclear Physics, Moscow State University, 119899, Moscow Russia, 
( 15 ) Nagoya University, Chikusa-ku, Nagoya, 4^4 Japan, 
( 16 ) University of Oregon, Eugene, Oregon 97403, 
Oxford University, Oxford, 0X1 3RH, United Kingdom, 
( 18 ^INFN Sezione di Perugia and Universita di Perugia, 1-06100 Perugia, Italy, 
^Rutherford Appleton Laboratory, Chilton, Didcot, Oxon 0X11 OQX United Kingdom, 

^Rutgers University, Piscataway, New Jersey 08855, 
( 21 ) Stanford Linear Accelerator Center, Stanford University, Stanford, California 94309, 
( 22 ) Soongsil University, Seoul, Korea 156-743, 
(23) xj n { vers Hy f Tennessee, Knoxville, Tennessee 37996, 
( 24 ) Tohoku University, Sendai 980, Japan, 
^University of California at Santa Barbara, Santa Barbara, California 93106, 
( 26 ) University of California at Santa Cruz, Santa Cruz, California 95064, 
( 2? ) Vanderbilt University, Nashville, Tennessee 37235, 
^University of Washington, Seattle, Washington 98105, 
( 29 ) University of Wisconsin, Madison, Wisconsin 53706, 
( 3 °) Yale University, New Haven, Connecticut 06511. 



(Dated: February 7, 2008) 



2 



Abstract 

A search for B® — B® oscillations is performed using a sample of 400,000 hadronic Z° decays 
collected by the SLD experiment. The B® candidates are reconstructed in the B® — > D~X channel 
with D~ — > (j)iT~ , K*°K~. The B® production flavor is determined using the large forward- 
backward asymmetry of polarized Z° — > 66 decays and charge information in the hemisphere 
opposite that of the B® candidate. The decay flavor is tagged by the charge of the Df. From a 
sample of 361 candidates with an average B® purity of 40%, we exclude the following values of the 
oscillation frequency: Am s < 1.4 ps _1 and 2.4 < Am s < 5.3 ps~ x at the 95% confidence level. 

PACS numbers: 13.20.He, 13.25.Hw, 14.40.Nd 
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I. INTRODUCTION 



The primary motivation for studying neutral B meson oscillations is to measure the 
poorly known Cabibbo-Kobayashi-Maskawa (CKM) matrix element V t d- The B® oscillation 
frequency corresponds to the mass difference, Am^, between the physical eigenstates of the 
B® — Bj system, which is sensitive to \Vtd\- Although Am^ is measured to within 2.5% [|l|, 
theoretical uncertainties lead to a 20% uncertainty in the extraction of \V t d\ [0- However, 
many uncertainties cancel in the ratio of mass differences in the B® and B® systems: 
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where and mg s are the B meson masses, fs d and /b s are the decay constants, and 
Bs d and B Ba are the "B-parameters" . Using this formula, and assuming |Vt s | = |Vd,|, one can 
obtain a 5-10% theoretical uncertainty on \V t d\ 0, ||- 

As yet, the B° s oscillation frequency has not been measured. The published lower limit on 
Am s based on the combined results from ALEPH, DELPHI, OPAL and CDF is 13.1 ps _1 
at the 95% confidence level |lj]. In the context of the Standard Model, other measurements 
suggest that Am s may be just beyond this current limit ||. 

This letter describes a study || of B° s — B® oscillations with the SLD experiment at 
SLAC. The measurement of B° s mixing requires excellent decay time resolution in order to 
resolve the very fast oscillations. The technique described in this letter, using B® — > D~X 
decays with D~ — > (f)7i~ or K*°K~ , has excellent decay length resolution (the best of any 
B° s analysis to date) and high purity (reconstruction of a D~ helps to discriminate against 
the prominent backgrounds from B + and B^ mesons). This work also represents the first 
use of polarized beams in a search for B° s mixing, which provide an effective new means of 
identifying the flavor of the B° s meson at production. 



II. APPARATUS AND EVENT SELECTION 



This analysis is based on a data set of 400,000 events of the form e + e~ — > Z Q — > hadrons, 
collected from 1996 through 1998 by the SLD experiment. A detailed description of the 
experiment can be found elsewhere [^|, |]. The analysis uses charged tracks reconstructed 
in the Central Drift Chamber (CDC) and the pixel-based CCD vertex detector (VXD3). 
The momentum resolution from the combined CDC and VXD3 fit is determined to be 
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c p ±/p± = 0.010 © 0.0024p_i_, where p± (in GeV/c) is the momentum of the track transverse 
to the beamline. The track impact parameter resolution in the transverse plane is a xy = 
8 © 33/(psin 3//2 9) \im and along the beam direction is a z = 10 © 33/(psin 3//2 9) fim, where 
p (in GeV/c) is the track momentum and 9 is the polar angle with respect to the electron 
beam. The tracking system is surrounded by the Cherenkov Ring Imaging Detector (CRID), 
a two-radiator system that allows good pion and kaon separation in the momentum range 
between 0.3 and 35 GeV/c. 

The goal of the analysis is to measure or constrain the oscillation frequency. This is 
accomplished by using B® decay candidates that are flavor tagged at both production and 
decay, and by measuring the proper time of the decay through measurements of the B° s decay 
length and energy. The decay time distributions of B mesons whose flavor at production 
and decay are different (the same), so called mixed (unmixed) events, are modulated by the 
oscillation frequency. A fit to these distributions constrains Am s . 

Event reconstruction consists of four main steps: hadronic event selection, Z° — > bb event 
selection, D~ decay reconstruction, and partial reconstruction of B® decays. A hadronic 
event is identified as having at least seven charged tracks, a total energy of at least 18 
GeV, and an event thrust axis satisfying \cos9 t h rust \ < 0.85. The thrust axis is calculated 
based on the energy clusters found in the liquid-argon calorimeter. The hadronic event 
selection removes essentially all dilepton events and other non-hadronic backgrounds. To 
enhance the fraction of Z° — > bb events in the sample, events are required to have at least 
one topologically reconstructed secondary vertex || with a vertex mass greater than 2 GeV 
in either hemisphere. The vertex mass calculation includes a correction for the missing 
momentum transverse to the B flight direction in order to partially account for missing 
particles. A neural network [[K|] is used to select a candidate secondary B vertex (if multiple 
topological secondary vertices exist) and its decay tracks. The resulting event sample is 97% 
pure bb with a single hemisphere b tagging efficiency of 54%. 

The D~ is reconstructed in one of two modes, D~ — > 07r~ or K*°K~ — > K + K~tt~. 
Oppositely charged tracks are first paired to form a <fi (K*°) candidate and a third track 
is then attached to form a D~ candidate. To maximize the discrimination between true 
D~ and combinatorial background events, kinematic information for the D~ candidate is 
fed into another neural network. In this case, the neural net inputs include CRID particle 
identification information for each of the three daughter tracks, the K + K~ (K~tt + ) invariant 
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FIG. 1: Distribution of KKir invariant mass for (fin and K*°K candidates combined, for the events 
that include a B® candidate vertex. The solid line shows a fit with the same functional form as is 
used in the likelihood function: two Gaussian distributions with common mean for the D s signal, 
a single Gaussian for the D + contribution and a second-order polynomial for the background. 

mass for cfi (K*°) candidates, the fit probability for the D~ decay vertex, the D~ decay length 
normalized by the decay length error, the total momentum of the Dj, the angle between 
the neutral meson (<fi or K*°) momentum in the D~ rest frame and the D~ flight direction, 
and the angle between the n~ or K~ from the D~ decay and the K + from the neutral 
meson decay in the rest frame of the neutral meson. The neural net cut that maximizes 
the sensitivity of the analysis to mixing is determined separately for each of the two D~ 
decay modes. The resulting KKn mass distribution for both modes combined is shown in 
Figure |l|. 

The B° s decay vertex is found by vertexing the virtual D~ track with other tracks in the 
hemisphere. The virtual D~ track is constructed by combining the 4-momenta of the three 
daughter tracks and constraining the parent D~ to pass through the decay vertex. The 
parent track error matrix is determined from propagation of the track measurement errors 
of the daughters. The B° s vertex fit is accomplished in two steps: 1) identify an intersection 
of the virtual D~ track with another charged track in the same hemisphere to act as a seed 
for the B® decay vertex and 2) add other charged tracks to the seed vertex if they are more 
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consistent with coming from the B s decay than coming from the primary interaction point 
(IP). To find the seed, the D~ track is individually vertexed with each track (excluding Dj 
daughters) in the same hemisphere. The vertex that is farthest from the IP and upstream 
of the D~ (or consistent with being upstream within 5cr) and has a vertex fit x 2 °f l ess 
than 5 is chosen as the seed. In order to determine if another track should be added to the 
vertex, we examine two parameters: 1) the distance D from the IP to the seed vertex and 
2) the distance L along the same direction from the IP to the point of closest approach of 
the track. If the ratio L/D is greater than 0.5 and the track forms a good vertex with the 
D~ (fit x 2 < 5), the track is added to the vertex. The latter condition is imposed to reject 
spurious tracks (often from a second charm decay if the B decayed to two charm particles) 
that do not point back to the B vertex. Only B decays for which the total charge of all 
associated tracks is Q=0 or ±1 are kept for the </>7r mode and those for which Q=0 are kept 
for the K*°K~ mode. The selected tracks are then vertexed together with the D~ to obtain 
the best estimate of the B decay position. 

The IP location, needed to calculate the B° s decay length, is determined from tracks 
that have vertex detector hits and extrapolate within 3a of the beamline. The coordinates 
transverse to the beam are averaged over approximately 30 sequential hadronic Z° events, 
while the position along the beam is determined event-by-event. The resulting error in IP 
position is 3.5/oti in the transverse plane and 17/xm along the beam, the best resolution of 
any high energy physics collider. By making use of this well-determined beam spot, the very 
small beam size, and information from the high precision CCD vertex detector, this analysis 
has a unique sensitivity for measuring decay times of the B mesons. 

An estimate of the B° s decay length resolution is determined event-by-event according to 
the vertex fit and IP uncertainties, with correction factors applied for each particle decay 
hypothesis (B®, B%, B + or 6-baryon) as determined from Monte Carlo simulation. The 
SLD Monte Carlo uses JETSET 7.4 with the B decay model tuned to CLEO and ARGUS 



data [11]]. Parameterizing the decay length resolution as a sum of two Gaussians with 
normalizations of 0.6 (core) and 0.4 (tail), the estimated multiplicative correction factors 
for the B® signal events are 1.07 and 2.16 for the core and tail resolutions, respectively. The 
resulting average decay length resolution for the B® signal events is 50 \im (core) and 151 
fim (tail). 

The B meson boost is calculated from separate estimates for the charged and neutral 
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particle contributions to the total energy. The charged energy is determined by summing 
all the charged tracks associated with the B decay assuming the pion mass (except for the 
two kaons from the D~ decay). The neutral energy estimate uses five different techniques. 
The first four techniques are calorimeter-based and use various constraints (beam energy, 
jet energy, B° s mass and calorimeter information) to estimate the neutral energy of the B 
|T2"|| . The fifth technique is based only on the information from the charged decay 



meson 



tracks and the kinematics of the decay (B vertex axis, charged track momentum and B® mass 
constraint) 0]. The results from the five algorithms are then averaged, taking correlations 
into account, to obtain the total B energy. The resulting average boost resolution (^f-) is 
represented by a sum of two Gaussians with widths (normalizations) of 8% and 18% (0.6 
and 0.4) for the B® signal events. 

Our final data sample includes 361 events within ±40MeV of the D~ mass peak (Figure [l]) 
with an average D~ purity of 48.1%. The composition of the D~ signal sample is calculated 
from the published branching ratio measurements with the relative reconstruction efficiencies 
of various decay modes taken from the Monte Carlo simulation. We estimate that the Dj 
signal peak consists of B° s -> DfX (55.1%), B° d -> DfX (22.4%), B + -> DfX (15.6%), 
6-baryon— > DfX (5.5%) and prompt c — ► D~ (1.4%). For the hadronic B° s decays, roughly 
10% of the decays yield a wrong-sign D s (B® — > DfX instead of 5° — > D~X), due to W + — > 
Df. Of the 361 B° s candidates, 39 are semileptonic decay candidates (B° s — > D~l + vX). The 
fraction of wrong-sign D s decays for the semileptonic modes is about 5%. These events 
have significantly higher B s fraction and tagging purity than the hadronic decay sample, 
and therefore are parameterized separately. 



III. TAGGING 

The flavor of the B® at decay is determined by the charge of the D~. A D~ is assumed 
to come from a B® and a Df is assumed to come from a B®. The B° s flavor at production is 
obtained by exploiting the large forward-backward asymmetry in polarized e + e^ R — ► Z° — > 
bb decays. The differential cross section for the decay is given by 

oc (1 - A e P e ){\ + cos 2 9) + 2A b (A e - P e ) cos9, (2) 
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where the asymmetry parameter can be expressed in terms of vector and axial- vector cou- 
plings, Af = 2afVf/(aj + vj) and 9 is defined as the angle between the outgoing fermion 
and the electron direction. The electron polarization is defined as P e = > where 

N(R) (N(L)) is the number of right-handed (left-handed) electrons in a beam bunch. The 
outgoing 6-quark is produced preferentially along the direction opposite to the spin of the Z° 
boson. Therefore, by knowing the polarization of the electron beam and the direction of the 
jet, the flavor of the primary quark in the jet can be inferred. The correct tag probability 
depends on the polar angle of the jet with respect to the incident electron beam direction and 
the electron beam polarization. The average electron beam polarization achieved during the 
run is about 73%. The resulting average correct tag probability is about 72%. In addition 
to the polarization tag, information from the hemisphere opposite to the reconstructed B° s 
is used to improve the identification of the production flavor. A series of neural networks 
is used to combine momentum- weighted jet charge, vertex charge, dipole charge ||, lepton 
charge and kaon charge information. The purity of the opposite hemisphere charged tags 
has been calibrated from the data to be (70.4±0.8)%. Combining all available tags, the 
overall production flavor tag purity is (77.8±0.8)%. In the end, candidates for which the 
flavor of the B meson at decay is more than 50% likely to be different from the flavor at 
production are said to be 'mixed'. 



IV. FITTING AND RESULTS 

To determine the B° s oscillation frequency, an unbinned likelihood function is used to 
describe the proper time distribution of mixed and unmixed events. For signal B° s — > DfX 
events tagged as mixed (unmixed), the proper time distribution has the form, 

P Bs (Am s ,t rec )^f xed = [ +00 -LeM-t/r Bs ) [ 1 =F (1 - 2tfc,J cos(Aro.f)] • 

e Bs (t)-G Bs (trec,t)-dt } (3) 

where t rec is the reconstructed proper time, t is the true proper time, r Bs is the B® lifetime, 
e Bg (t) is the reconstruction efficiency and G Bg (t rec ,t) is the proper time resolution function 
for the B° s events. The overall mistag probability is r] Bg = r]i(l — r] B J + (1 — rji)r] Bs , where r\i 
is the production flavor mistag probability and rj Bg is the decay flavor mistag probability for 
the B® events. The efficiency and proper time resolution functions are derived from Monte 
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Carlo simulations. The Monte Carlo vertexing resolution is based on our understanding of 
the impact parameter resolution, which is carefully tuned to match the data. A comparison 
of decay length resolution between data and Monte Carlo using 3-prong r decays shows 
good agreement. The proper time resolution is expressed in terms of decay length and boost 
resolutions (o"L,cr 7/ g) according to: 

1/2 



j(3c \ 7/3 



(4) 



where the indices i, j =1 for core and 2 for tail resolution. Both <jl and a 1 p are determined 
event-by-event. The decay length resolution <jl depends on the B° s vertex fit error matrix and 
IP uncertainties, and the boost resolution a 1 p is parameterized as a function of charged track 
energy in the decay. The decay time resolution function is comprised of four at Gaussian 
distributions given by the various and cr 7j g core-tail combinations. 

The distribution for B® mesons, which are also subject to oscillations, is identical to 
equation |3] with B° s subscripts replaced by B® subscripts. The terms for B + and 6-baryon, 
which do not oscillate, are constructed in an analogous fashion, but with Am = 0. The 
probability density function for the data sample has the form, 

+ [1 — /dJ Pcomb- (5) 

The fraction of D~ signal above the combinatorial ttikkit background, fp s , is estimated from 
the previous D~ mass fit (Figure |IJ) as a function of rriKKn- The remaining fractions are 
calculated based on the measured branching ratios with reconstruction efficiencies taken from 
the Monte Carlo simulation. The function P cc describes the prompt charm (Z° — > cc) events, 
and is obtained from Monte Carlo simulation. The time distribution for the combinatorial 
events, P com b is parameterized directly from the data using events in the D~ mass sidebands. 
The sideband regions are defined as: 1.7 < rriKKn < 1.8 GeV (lower sideband) and 2.05 
< rriKK-K < 2.2 GeV (upper sideband). The parameters iVj, are normalization constants, 
obtained by integrating the sum of Pmixed and P un mixed over all reconstructed proper time. 



In the absence of a signal, the amplitude fit method |14[ is used to set a limit on Am s . 
The amplitude fit is equivalent to a Fourier analysis, in which one searches for peaks in the 
frequency spectrum of oscillation. To perform an amplitude fit, the likelihood function is 
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FIG. 2: Measured B® oscillation amplitude as a function of Am s . The light-grey (dark- grey) band 
shows the 90% confidence region obtained from statistical (total) uncertainties. Values of Am s for 
which the allowed band is below A = 1.0 are excluded at the one-sided 95% confidence level. 



modified by replacing cos(Am s t) with A ■ cos(Am s t). The amplitude A and its error a a 
are then measured at each assumed value of Am s . If mixing occurs at the chosen value of 
Am s , the fitted value of A should be consistent with unity. At values of Am s sufficiently 
far from the true mixing frequency, the fitted value of A should be close to zero, consistent 
with no oscillation. Values of Am s for which A + 1.645cta < 1 can be excluded at the 95% 
confidence level. We have tested this procedure on simulated data for Am s values of 4, 10, 
17 and 270 ps _1 to verify that the amplitude fit behaves as expected in all cases. 

The amplitude plot for this analysis is shown in Figure ||. There is no evidence for a 
significant signal anywhere in the plot, so we use the analysis to set a limit on Am s . The 
systematic errors on A are evaluated according to reference ||14j| , all of which are negligible 
compared to the statistical errors. The systematics are dominated by uncertainties in a small 
reconstructed proper time bias (evaluated at 100% of the simulated correction, typically a few 
hundreths of a picosecond), in the production flavor tag (evaluated with ±0.8% uncertainty 
in the average tag probability and a reweighting of the shape of the distribution as a function 
of neural net output), in f D - (as determined from the KK-k mass fit), and in the boost 
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TABLE I: B production fractions and various branching ratios assumed in the amplitude fit. The 
uncertainties for the branching ratios do not include the uncertainty from B(D S — > <fiir). 



Parameter 



Value and Error 



Ref. 



f(b - B° a ) 
f(b^B%B+) 
f(b — > b — baryon) 



0.100 ±0.012 
0.401 ±0.010 
0.099 ±0.017 















R b .f(b^W s )- B(B° -> £>+X) • B(D+ - 
/(6 _^ W - _ D -) . B(D- -> 07T-) 
B{B^,B+^DfX)-B{D-^<fm-) 

D-X)/B(B° d ,B+^DfX) 



■7T 



(6.211°;?*) x 10~ 4 
(3.66 ± 0.45) x 10- 3 
(3.71 ± 0.28) x 10-3 
0.172 ±0.083 

(3.4 ±0.3) x 10~ 3 



0] 
0] 



0] 



0] 



resolution (10%). To a lesser extent, there are contributions from uncertainties in decay 
length resolution (7%) and branching fractions assumed in the fit. Uncertainties in particle 
lifetimes and in the B® oscillation frequency are completely negligible. The list of systematic 
errors is shown in Table [TT]. The dark-grey band in Figure 0, barely visible at the edge of the 
light-grey band, shows the effect of adding the total systematic error to the calculations. 

Based on the result of the amplitude fit, the values of the B° s oscillation frequency excluded 
at the 95% confidence level are: Am s < 1.4 ps" 1 and 2.4 < Am s < 5.3 ps -1 . 

V. CONCLUSIONS 

In conclusion, we have performed a search for B° s — B® oscillations using 361 B° s — > D~X 
candidate events with an average B° s purity of 40%. The excluded values of oscillation 
frequency are: Am s < 1.4 ps" 1 and 2.4 < Am s < 5.3 ps -1 at the 95% confidence level. 
The analysis exploits the unique large polarized forward-backward asymmetry in Z° — > bb 
decays to enhance the production flavor tag. Combining the good production flavor tag and 
excellent proper time resolution, the analysis contributes to the world average at high Am s 
despite low statistics. Taken in conjunction with other published results, our result raises 
the 95% C.L. world limit on Am s from 13.1 ps -1 to 13.9 ps -1 . 
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TABLE II: Table of statistical and dominant systematic uncertainties for several Am s values. 



Am s 


10 ps -1 


15 ps" 1 


20 ps" 1 


Measured amplitude A 


0.029 


1.027 


2.513 


—stat 
°A 


±0.933 


±1.361 


±2.283 


syst 
°A 


+0.088 
-0.084 


+0.312 
-0.313 


+0.776 
-0.789 




+0.012 
-0.012 


+0.021 
-0.022 


+0.010 
-0.011 


f(b — > b— baryon) 


+0.004 

-o!oo4 


+0.001 

-o!ooi 


—0.019 

+o!oi9 


R b -B(b^B°)-B(B°^D+X) 


— U.Uoz 
+0.043 


n ha n 
— U.U4U 

+0.048 


— U.UDo 
+0.068 


B(b^W-^D-) 


+0.013 
-0.013 


+0.022 
-0.022 


+0.010 
-0.012 


B(B%B+^DfX) 


+0.013 
-0.014 


+0.013 
-0.014 


+0.042 
-0.047 


B(B°, B+^D7X) 
B{B%B+^DfX) 


-0.013 


-|-U.UOo 
-0.060 


-0.156 


B(c^D-) 


+0.003 
-0.003 


+0.007 
-0.007 


+0.011 
-0.011 


Decay length resolution 


+0.034 
-0.036 


+0.015 
-0.015 


+0.026 
-0.045 


Boost resolution 


+0.049 
-0.048 


-0.052 
+0.042 


-0.161 
+0.140 


Jd s 


+0.019 
-0.018 


+0.096 
-0.096 


-0.090 
+0.076 


Average production flavor tag 


-0.029 
+0.031 


-0.036 
+0.036 


-0.051 
+0.048 


Production flavor tag shape 


-0.020 
+0.018 


+0.008 
-0.013 


+0.306 
-0.318 


Proper time offset 


+0.007 
-0.007 


+0.278 
-0.278 


+0.671 
-0.671 
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